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Abstract 

We have analyzed the quantum interference effects in the macro- 
scopic "superconducting molecule". The composite system consists of 
two massive superconducting rings, each interrupted by a Josephson 
junction, which are at the same time weakly coupled with one another. 
The special case of coupling via the Josephson four-terminal junction 
is considered. The structure of the macroscopic quantum states in an 
applied magnetic field is calculated. It is shown, that depending on the 
values of the magnetic fluxes through each ring, the system displays 
two groups of states, the " orthostates" with both induced currents 
going in the same direction, and the "parastates" with the opposite 
currents and with the total induced flux locked to zero value. The 
transition to the flux locked state with changing of the total applied 
flux is sudden and is preserved in a certain interval which is determined 
by the difference of the fluxes applied through each ring. It makes the 
system sensitive to small gradients of the external magnetic field. 
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The system which we studied, is shown in Fig.l and consists of two bulk 
superconducting rings, coupled via the 4-terminal Josephson junction §]. 
The 4-terminal Josephson junction is a system of two microbridges, 1 — 2 and 
3 — 4, having the common centre "o". The interference in the cross section 
"o" of macroscopic wave functions of the jth terminal (j = 1, ..4) leads 
to nonlinear coupling and consequently to interference between the current 
states in each ring. The resulting current state of the whole system can 
be regulated by the difference of the magnetic fluxes applied through the 
rings, in analogy with the phase difference between two weakly coupled bulk 
superconductors. The studying of the macroscopic quantum states of such 
"superconducting molecule" is the subject of the present paper. 

The free energy U of our system in an applied magnetic field contains the 
magnetic energy U m and the Josephson coupling energy Uj. The energy U m 
has the form ||: 



U m = l> 2 + ±12 2 > 1 ($|-$ 2 ), (1) 

1\L\Li — L l2 ) 2{L\L 2 — L l2 ) {L\L 2 — L l2 ) 

where §\ 2 are the external magnetic fluxes applied to the rings 1, 2 and <3>i ; 2 
are the resulting fluxes embraced by the rings; and L\ 2 are the ring self- 
inductances and the mutual inductance [L\ 2 < LiL 2 ) . The coupling energy 
Uj (in dimensionless units) is expressed in terms of phases ipj (j = 1, ..4) of 
the superconducting order parameter in the jth terminal : 

rr 2 2^1 2^2 „ <Pl <p2 (r) s 

U j = —k cos cos 2k cos — cos — COS Y, z 

2 2 2 2 w 

if we introduce the phase differences across the weak links in the rings 

01 = <fi - ¥2, 4>2=<-P3- V?4 

and the "total" phase difference between the rings 

_ <fl + <f2 _ ^3 + V?4 

X ~ 2 2 
The coupling constant k in (2) is the ratio of critical currents of the weak 
links 1 — 2 and 3 — 4. In the following for simplicity we will consider the case 
of identical rings with L\ = L 2 = L and the symmetrical coupling k — 1 

Uc,12 — -^c,34 = Ic)- 
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The phase differences 0i 5 2 are related to the magnetic fluxes $1,2 by : 
</?i ; 2 = — 2e$i j2 /^. Thus, the total energy in reduced units of the two coupled 
rings as function of the embraced magnetic fluxes at given values of the 
applied fluxes is defined as 



$ 2 , x 1 *!, *5) = ^ + -m- - **)- 

(3) 

— cos 2 4^ — cos 2 ^ — 2 cos ^ cos ^ cos x, 

where £ = L± 2 /L the normalized mutual inductance (£ < 1), £ = {2eI c /%)L{l— 
£ 2 ) the dimensionless effective self-inductance; the magnetic fluxes are mea- 
sured in units h/2e. Note the dependence of the potential U on the phase x- 
As we will see, in the stable steady state the phase x can take only the value 
or 7r, which corresponds to existence of two groups of states with different 
symmetry. 

The minima of the potential U (3) with respect to variables $1, $2) X 
at given external fluxes and $2 determine the stable steady states of our 
system. The minimization of U with respect to x gives that the phase x 
takes the value or n, depending on the equilibrium values of $iand $2 

cos% = sign(cos cos -^-). (4) 

In the steady state = -j^. = 0, or : 

$1 $1 $2 
$1 - £$2 = $1 - + £ sin— [cos — + cos x cos — ] , (5) 



-£<&l = $ 2 + £ sin — [cos — + cosxcos—)], (6) 

2 2 2 

with cosx defined by the condition (4). 

The solutions of eqs.(5) and (6) {$1, $2} which correspond to the minima 
of the potential U must satisfy the requirements: 

d 2 U d 2 U d 2 Ud 2 U , d 2 U 



d&( > °' a$i > °' d$ 2 &s>i >0 ' (7) 

It can be shown that the conditions (7) are fulfilled for all values of <£>iand $2 
if £+£ < 1. In the following we will consider the case when the inductances C 
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and I are small enough to satisfy this inequality. Thus, all solutions {$1, $2} 
of the equations (4), (5) and (6) determine the possible stable or metastable 
states of the system. The circulating ringcurrents Ii j2 in state {$1, $ 2 } are: 

1 $1 $1 $2 

I\ — — sin $! — sin —sign(cos — ) | cos — |, (8) 
2 
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J 2 = — - sin $ 2 — sin -^-sign(cos —^) | cos | (9) 
in units of J c . 

The value of cosx in eqs.(5), (6), which equals ±1, determines two pos- 
sible "binding" of the current states in individual rings. The group of states 
{$i,$2,X = 0} we will call symmetric, or "ortho", states and the group of 
states {$1, $ 2 , x — - antisymmetric, or "para", states. As we will see, the 
first one corresponds to the induced ringcurrents going in the same direction, 
and the second one - to the currents going opposite. 

We will study the behaviour of our system in an applied magnetic field 
as the response on the total applied magnetic flux $ e = $^ + $2 a ^ given 
difference 5 e = <&\ — $| of the fluxes through the each ring. The state of 
the system as whole is determined by the total embraced magnetic flux $ = 
$1 + $2 , or by the total orbital magnetic moment M, which is proportional 
to the sum of the induced ringcurrents , I = I\ + 7 2 . Note, that the positive 
(negative) sign of I corresponds to the parallel (antiparallel) direction of M 
with respect to the external magnetic field H. From the (4-6) we obtain: 

C $ 5 
$e = $+ l37 sin 2"f COS 2 +C ° S ^' ^ 

£ 5 $ 
5 e = 5 + Y^sin-[cos- + cosx], (11) 

$ 5 

cos x = sign(cos — + cos -), (12) 

where 5= $1 — $2- 

Let us start from the case of small inductances £, C <C 1. If 5 e = 0, from 
the eqs.(ll), (12) follows that 5 = and x = 0. For $($ e ) we have the 
usual equation $ e = $ + 2£sin^ for the case of decoupled rings j4j, each 
interrupted by a Josephson junction. At 5 e 7^ and consequently 5^0, the 
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solutions with % = 7r are possible in the vicinity of $ ~ 2tx. In the limit 
£ — > for the total induced magnetic flux $($ e , 5 e ) we have the expression 



$ = $ e — £ sin — [cos — + sign(cos — + cos — )]. 



(13) 



In the case of small 5 e <C 1 it becomes : 



$ e - 2£sin^ 


if 




- 2tt 


l>l ^ 1 




if 




- 2tt 





(14) 



Thus, for given value of S e with changing of the total applied flux $ e the 
system switches from the state with x = to the state with \ — n - I n 
interval 2n — 5 e < $ e < 2tt + 5 e the total induced flux $ — $ e equals to zero 
for 5 e <C 1. We call such behaviour magnetic flux locking. It is emphasized 
that the transition to the flux locked state is sudden and is preserved in a 
certain interval of the applied magnetic flux. For the sum of the induced 
ringcurrents I — I\ + 12, in the limit C — > we have 



In the flux locked state the total current I equals to zero in correspondence 
with (14). Thus the ringcurrents are going in opposite directions and 
compensate each other, or the system is in the "para" state . The complete 
compensation takes place for 5 e 1, with the corrections to zero value being 
of the order of (<5 e ) 2 . In Fig. 2 we plot the dependence of 7($ e ) (15) for the 
flux difference 5 e = 27r/10, or in dimension units 1/10 of a flux quantum 
h/2e. The dashed line is the sum of the currents in the same, but decoupled, 
rings with the same applied fluxes <&\ = l/2($ e + 5 e ) and $2 = l/2($ e - 5 e ) 
. The magnetic susceptibility of the system as whole is proportional to — 
and will reflect the behaviour of the induced currents. 

For finite, but small, values of the inductances , the behaviour described 
above will be qualitatively the same. Only instead of the sharp switches 
hysteretic regions appear, of which the width is proportional to £ . In Fig. 3 
the dependence $($ e ) for £ = 0.25, £ = and 5 e = 1 is shown, as follows from 
the numerical solution of the eqs. (10-12). Naturally, these hysteretic regions 
will be smeared by thermal fluctuations (see the analysis of the influence of 
noise on the similar system, the so called 4-terminal SQUID, in ref. 0). 



J($ e ,5 e ) 



sin — [cos h sign(cos h cos — )]. 



(15) 
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In conclusion, we have studied the macroscopic quantum states in the 
system of two weakly coupled superconducting rings. The nonlinear coupling 
leads to interference between the current states in each ring. It is manifested 
as the cooperative behaviour of the rings in some region of the applied mag- 
netic fluxes, which we call magnetic flux locking. We would like to remark 
that our macroscopic approach is not restricted by the special kind of the 
coupling through the crossed superconducting bridges. In fact, any meso- 
scopic 4-terminal weak link will produce a coupling similar to the Uj (2). 
For example, it can be the experimental setup described in ref.|§, namely 
the two-dimensional normal layer which is connected with four terminals 
instead of the two ones as studied in ref.||. 

One of the authors, A.N.O., would like to acknowledge the support for 
this research from the Kamerlingh Onnes Laboratory, Leiden University. 
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FIGURE CAPTIONS 

Figure 1. The two bulk superconducting rings, coupled via the 4- 
terminal Josephson junction ( the region closed by the dashed lines, of which 
the area is of the order of the coherence length squared). 

Figure 2. The total induced current as a function of the total applied 
flux at given difference of applied fluxes through each ring 5 e = 27r/10 (or 
1/10 of a flux quantum h/2e). C = 0. The dashed line is the corresponding 
dependence in the case of decoupled rings. 

Figure 3. The dependence of the total magnetic flux $ on the total 
external flux $ e for 5 e = 1, C — 0.25, £ — 0. The arrows indicate the jumps 
of the flux from metastable to stable states. The dashed line is $ = <3> e . 
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